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ABSTRACT 

In this work we present IRAM-30m telescope observations of a sample of bulge- 
dominated galaxies with large dust lanes, which have had a recent minor merger. We find 
these galaxies are very gas rich, with H 2 masses between 4xl0 8 and 2xl0 10 M 0 . We use 
these molecular gas masses, combined with atomic gas masses from an accompanying paper, 
to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample 
objects vary widely (between «50 and 750), suggesting many objects have low gas-phase 
metallicities, and thus that the gas has been accreted through a recent merger with a lower 
mass companion. We calculate the implied minor companion masses and gas fractions, find¬ 
ing a median predicted stellar mass ratio of «40:1. The minor companion likely had masses 
between ^10 7 - 10 10 M 0 . The implied merger mass ratios are consistent with the expectation 
for low redshift gas-rich mergers from simulations. We then go on to present evidence that 
(no matter which star-formation rate indicator is used) our sample objects have very low star- 
formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type 
galaxies from ATLAS 3D which already show a suppression. This suggests that minor mergers 
can actually suppress star-formation activity. We discuss mechanisms that could cause such a 
suppression, include dynamical effects induced by the minor merger. 

Key words: galaxies: elliptical and lenticular, cD - ISM: molecules - galaxies: ISM - galax¬ 
ies: evolution - galaxies: interactions 


1 Introduction 


Early type galaxies (ETGs) have traditionally been thought of 
as gas-poor objects. This perception has slowly been challenged 
(e.g.|Roberts et al.|1991|), and recent results (e.g. Combes, Young & 


Bure au|2007||Welch, Sage & Young|2010| | Young et al.|2011| > have 
showed with statistically complete samples that ^1/4 of ETGs have 


* E-mail: t. davis4 @herts . ac . uk 


>10 7 M 0 of cold molecular gas, and ^40% of ETGs have sizeable 
atomic gas reservoirs ( [Morganti et al.|20 06 , Sage & Welch| 2006 , 
Serra et al.|2012| >. Dust is also present in a large fraction of ETGs, 
e.g. |Smith et al.f < [2~012| ) find that ^50% of ETGs in the Herschel 
Reference Survey have cold dust. Low level residual star forma¬ 
tion has also been detected through studies of UV emission (e.g. 
Yi et al.|2005||Kaviraj et al.|2007||Salim & Rich|2010||Wei et al.| 

2010| ), optical emission lines (e.g. |Crocker et al.|2011|) and infrared 
emission (e.g. |Knapp et al.| 1989) |Combes, Young & Bureau|2007| 
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Temi, Brighenti & Mathews 2009, Shapiro et al. 2010, Davis et al. 


20l4j. 


Observationally, in field environments, the gas and dust in 
ETGs seems to come mainly from external sources. The evidence 
for this comes from the kinematic decoupling of the cold and ion- 
ized gas from the stars (|Sarzi et al.|2006}|Davis et al.|201 l~a||Katkov^] 
|SiTchenko & A fanasiev 2 014| ) , the orders of magnitude surplus of 
ISM compared to what one would expect from stellar mass loss (e.g. 


Merl uzzi|1998||Knapp et al.|T 989 Row lands et al.|201~2l|Kaviraj1 


et al.|2012), he large spread in gas-to-dust ratios within their cold 
ISM ( [Smith et al.||2012| ), the large fraction of gas rich ETGs that 
are morphologically disturbed in deeper imaging (e.g. jvan Dokkum| 
|2005||Puc et al.|20T5] > and the presence of young kinematically de¬ 
coupled cores in IFU surveys (e.g. |McDermid et al.|2006| ). 

The external sources that supply this gas likely include both 
mergers with gas-rich objects and accretion of gas from the inter- 
galactic medium. Major mergers are thought to be too rare in the 
low redshift universe to supply enough gas, so minor mergers are 
argued to be the dominant source of gas (e.g. [Kaviraj et al.|2009] 
|2011| ). Indeed, some studies claim that minor mergers are likely 
to be driving at least half of the cosmic star formation budget in the 
local Universe ( [Kaviraj |2014a|b| ). We thus need a better understand¬ 
ing of the role of minor mergers, in ETGs in particular, and galaxy 
evolution in general. 

The nearest large ETG to us, Centaurus A, is a good exam¬ 
ple of this process in action. Its 10 9 M© of cold gas and 10 7 M© of 
dust were likely acquired in a recent, relatively minor merger ( [Israeli 
1998; Parkin et al. 2012). This gas/dust reservoir makes its presence 
obvious in optical images as a large warped dust lane, seen in ab¬ 
sorption against the stellar continuum. By selecting disturbed ETGs 
in the field with such large dust lanes we are able to probe the minor 
merging process, which supplies the fuel that is rejuvenating these 
objects. 

Theoretically it has been argued that the main mode of gas 
supply in ETGs in the local universe is accretion of gas from the 
inter-galactic and circum-galactic medium ( [Lagos et al.|2014| ), with 
minor mergers playing a secondary, albeit still very relevant role. 
Understanding if the observational and theoretical models for gas 
accretion onto massive galaxies can be reconciled is an important 
area of future research (e.g. |Lagos et al.|20l5] >. 

Dust-lane ETGs are relatively rare and the properties of the 
gas in such objects have yet to be studied in a statistical sample. 
In this work we present molecular gas observations of a sample of 
bulge-dominated galaxies with large far infrared bright dust lanes, 
selected from the unbiased Sloan Digital Sky Survey (SDSS) pho¬ 
tometry in Kaviraj et al. (2012 ) as part of the Galaxy Zoo project 
( |Lintott et al.||2008] r See Section [2] for the full selection criteria. 
These objects are all recent minor merger candidates (Oosterloo 
|et al.||2002l |Shabala et al.||2012| ). We selected the most dust rich 
objects (using 250/xm luminosities) and observed them with the In- 
stitut de Radioastronomie Millimterique (IRAM) 30m telescope, in 
order to estimate the mass of molecular gas present in each system. 

In this paper we present the IRAM-30m observations (see Sec¬ 
tion [3j, and the derived quantities (see Section [4). In Section [5] we 
discuss the dust-to-gas and dust-to-stellar mass ratios of these sys¬ 
tems, derived from Hersche^ space telescope ( [Pilbratt et al.|2010| 
observations from the Herschel -ATLAS survey (H-ATLAS; |Eales| 


f H erschel is an ESA spac e obse rvatory wit h science i nstruments provided 
by European-led Principal Investigator consortia and with important partic¬ 
ipation from NASA. 


|et al.|2010) . We use these quantities to constrain the mass of the 
dwarf companion systems (and thus the mass ratio of the merg¬ 
ers). We also explore how the minor merging process affects the 
efficiency with which gas turns into stars. Finally we conclude in 
Section[6] In this paper we have assumed a cosmology with Ho= 71 
km s- 1 Mpc -1 , 0.27 and Ua=0.73. 


2 Sample 

The [Kaviraj et al.| ( |2012[ hereafter K12) dust lane sample is 
drawn from the entire SDSS Data Release 7 (which covers 11663 
deg 2 ) through visual inspection of galaxy images via the Galaxy 
Zoo (GZ; |Lintott et al.||[2008| ) project. Approximately 19,000 GZ 
galaxies at z < 0.1 were flagged as containing a dust lane by at 
least one GZ user (individual galaxies have 50+ classifications), and 
each galaxy was re-inspected by our team, as described in K12. 
This yields a final sample of 352 dust lane galaxies, with a median 
redshift, absolute r-band magnitude and stellar mass of 0.04, -21.5 
and 6xl0 10 M© respectively. These objects are massive, have red 
optical colours and are bulge-dominated. Due to the Galaxy Zoo 
selection criteria, classical ellipticals, lenticular galaxies, and some 
massive bulge-dominated Sa-type spirals are included. Despite this, 
we refer to these objects as dust lane early-type galaxies (DETGs) 
in this work to be consistent with previous studies. 

The DETGs in this study were extracted from this parent sam¬ 
ple by cross-matching the K12 objects with the Phase 1 and North 
Galactic Pole (NGP) fields of the Herschel -ATLAS survey (H- 
ATLAS; [Eales et ah||2010| ), as described in |Kaviraj et al.| ( [2013[ 
hereafter K13). The final DETG sample of K13 comprises 23 galax¬ 
ies, with a median redshift of 0.05 and stellar mass of 8xl0 10 M 0 
respectively. The sample of DETGs discussed in this paper are the 
17 dustiest objects (estimated using the 250-500/mi data) from the 
K13 sample. These objects are all recent minor merger candidates, 
having clear tidal low surface brightness features visible in their 
optical images (see |Shabala et al.|2012| for full details). 

2.1 Galaxy physical parameters 


Twenty filter photometry is available for each of these objects, cov¬ 
ering three orders of magnitude in wavelength (0.12 - 500 /xm;|Ibar| 


et al.|20lO[ Pascale et al.|2011||Rigby et al.|2011||Smith et al.|2011| 

Driv er et al.||201 1 These data are used in K13 to estimate the 

physical properties of our DETGs using the energy-balance tech¬ 
nique of [Da Cunha, Chariot & Elbazj < [2008[ DC08 hereafter) via 
the MAGPHYS cod^f] We refer the reader to K13 for more details, 
but outline below the basic premise of these models. 

In each model it is assumed that the UV-optical radiation emit¬ 
ted by stellar populations is partially absorbed by dust, and this ab¬ 
sorbed energy is matched to that re-radiated in the FIR. The optical 
library of 50,000 spectra is produced using the the latest version of 
the population synthesis code of Bruzual & Chariot (2003), Char¬ 
iot & Bruzual (2007, in prep), and assumes exponentially declining 
SFHs with superimposed stochastic bursts. The model spectra cover 
a wide range of age, metallicity, SFH and dust attenuation and as¬ 
sume a [Chabrier] ( |2003j ) IMF. The infrared libraries contain 50,000 
SEDs comprised of four different temperature dust components, 
from which the dust mass is calculated. In stellar birth clouds, these 
components are polycyclic aromatic hydrocarbons (PAHs), hot dust 
(stochastically heated small grains with a temperature 130 — 250 K), 
and warm dust in thermal equilibrium (30 — 60 K). In the diffuse 
ISM the relative fractions of these three dust components are fixed, 


J The DC08 models are publicly available as a user-friendly model package 
MAGPHYS atwww.iap.fr/magphys/. 
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Table 1 . Properties of the galaxies observed in this work 


NED name 

DETG ID 

z 

Dlum 

M* 

M d 

SFRmAGPHYS 




Mpc 

logio(M 0 ) 

logio(M 0 ) 

(Meyr- 1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

2MASXJ09033081-0106127 

18 

0.040 

174.2 

10.9d= 0.1 

7.41± 0.05 

0 09+ u * U5 

u - u ^-o.oi 

2MASXJ14161186+0152048 

106 

0.082 

368.3 

11.2d= 0.1 

8.11d= 0.03 

0 02 +0 ' 06 
U,U -0.01 

CGCG018-061 

231 

0.030 

129.7 

10.5± 0.2 

7.65± 0.08 

0 65 +0 ‘ 15 

U - OJ -0.17 

2MASXJ13341710+3455455 

920 

0.025 

107.6 

10.8±0.1 

7.28=b 0.07 

0 05 +o ‘° 3 

U ’ UJ -0.03 

NGC5233 

967 

0.026 

112.0 

ll.liO.l 

7.59± 0.08 

0.02±o ; °t 

2MFGC10709 

984 

0.037 

160.8 

10.8±0.1 

7.43=b 0.07 

°- 08 to.04 

UGC08505 

1004 

0.040 

174.2 

10.8±0.1 

8.12± 0.06 

O Q7+0.40 
-0.50 

UGC08536 

1014 

0.025 

107.6 

10.6± 0.1 

7.43=b 0.06 

1 OQ+0-41 
i,Z °-0.43 

2MASX J13194502+3134277 

1030 

0.045 

196.7 

10.7d= 0.1 

7.41± 0.03 

ooi±g;gi 

2MASXJ13411704+2616190 

1111 

0.064 

283.7 

ll.liO.l 

8.18± 0.07 

0-68^0 29 

NGC4797 

1232 

0.026 

112.0 

11.2± 0.1 

7.00± 0.06 

0 03 +0 ‘ 03 
U - U:> -0.01 

2MASXJ12581769+2708117 

1265 

0.040 

174.2 

10.8±0.1 

7.45=b 0.08 

0 09 +0 - 03 
U ’ Uy -0.02 

2MASXJ13010083+2701312 

1267 

0.078 

349.3 

10.9d= 0.1 

7.95± 0.08 

0 44+ 0 - 44 
U/f ^-0.37 

2MASXJ13333299+2616190 

1268 

0.037 

160.8 

10.8±0.1 

7.04 =b 0.03 

0 01+°- 06 
U ’ U1 -0.01 

2MASX J12561089+2332394 

1362 

0.133 

618.6 

11.0± 0.1 

8.24=b 0.03 

2 51+ 1 - 48 
z,ja -1.50 

2M ASX J13244786+2417386 

1372 

0.065 

288.4 

10.7d= 0.1 

7.80± 0.08 

0.231°;^ 

CGCG013-092 

1947 

0.035 

151.9 

10.9d= 0.1 

7.79=b 0.06 

O.28l° 0 : 0 04 


Notes: Column 1 shows the NED designator for these objects. Column 2 shows the DETG ID number for these sources from K13. Column 3 
and 4 show the spectroscopic redshift, and associated luminosity distance (assuming Ho= 71 km s -1 Mpc -1 , Q m = 0-27 and Q^=0.73). The 
stellar mass, dust mass and star-formation rate (averaged over the last 10 8 years) of these objects are in Column 5, 6 and 7, derived as in K13. 


but an additional cold dust component with an adjustable tempera¬ 
ture between 15 and 25 K is added. The dust mass absorption coef¬ 
ficient oc A _/3 has a normalisation of Kgso = 0.077 m 2 kg -1 . 
A dust emissivity index of (3 = 1.5 is assumed for warm dust, and 
P = 2.0 for the cold dust component. 

The physical properties of our sample DETGs are presented 
and discussed in detail in K13, and here we use similar values de¬ 
rived using the same modelling code, but with updated PACS and 
SPIRE photometry (Valiante et al., in prep). The galaxy properties 
(and common names extracted from the NASA/IPAC Extragalactic 
Database; NEE[§]) are presented in Table[l] We note that the changes 
caused by using the updated photometry are generally small, and 
the results discussed here are robust to the inclusion of this repro¬ 
cessed Herschel data. 

3 IRAM-30m Data 

The IRAM 30-m telescope at Pico Veleta, Spain, was used 
between the 4th and 7th of January 2013 to observe CO emis¬ 
sion in our sample of dust rich DETGs (proposal 182-12, PI 
Davis/Kaviraj). We aimed to detect CO(l-O) and CO(2-l) simul¬ 
taneously in the 3mm and 1mm atmospheric windows. The beam 
full width at half-maximum (FWHM) of the IRAM-30m at the fre¬ 
quency of these lines is 21'/3 and 10'/7, respectively. The Eight 
Mixer Receiver (EMIR) was used for observations in the wobbler 
switching mode, with reference positions offset by d=100" in az¬ 
imuth. The Fourier Transform Spectrograph (FTS) back-end gave 
an effective total bandwidth of ~4 GHz per window, and a raw 
spectral resolution of 200 kHz (~0.6 km s -1 at A =3mm, w0.3 
km s- 1 at A =1 mm). The Wideband Line Multiple Autocorrela¬ 
tor (WILMA) back-end was used simultaneously with the FTS, and 
WILMA data were used at 230 GHz when possible to avoid known 
issues with platforming in the FTS. 

The system temperatures ranged between 90 and 200 K at 
3 mm and between 100 and 220 K at 1 mm. The pointing was 

§ http://ned.ipac.caltech.edu/ 
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checked every 2 hours on a nearby planet or bright quasar, and 
the focus was checked at the beginning of each night, after sun¬ 
rise, or more often if a suitable planet was available. The time on 
source ranged from 15 to 190 min, being weather-dependent, and 
was interactively adjusted to try and ensure detection of molecular 
emission. 

The individual ^6 minute scans were inspected, and baselined, 
using a zeroth-, first- or second-order polynomial, depending on the 
scan. Scans with poor quality baselines or other problems were dis¬ 
carded. The good scans were averaged together, weighted by the 
inverse square of the system temperature. We consider emission 
lines where the integrated intensity has greater than a 3cr signifi¬ 
cance (including the baseline uncertainty; |Crocker et al. 2 012) to 
be detected. 

Integrated intensities for each detected line in each galaxy 
were computed by fitting a Gaussian profile in the class package 
of gilda^| or summing the spectrum over the velocity (where 
the profile shape is non-Gaussian). Both methods produce consis¬ 
tent results, and here we present the values derived from the sum. 
We convert the spectra from the observed antenna temperature (T*) 
into units of Jy, utilising the point source conversion S Jy / K as tab¬ 
ulated on the IRAM websit^jj] Where the exact frequency has no 
measured efficiency available we linearly interpolated the values 
required. Table [2] lists the RMS noise levels in each spectrum, the 
velocity width summed over, and the line integrated intensities (in 
Jy km s -1 ). Converting these integrated line intensities to the main 
beam temperature scale can be accomplished by dividing by a factor 
of 1.28 at 3mm§. 


% http://www.iram.fr/IRAMFR/GILDAS/ - accessed 14/01/13 

|| http://www.iram.es/IRAMES/mainWiki/EmirforAstronomers - accessed 

23/04/2013 
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Figure 1 . Three-colour (gri) optical images extracted from the SDSS (each with an identical angular size of 50x50"; left column), as well as CO(l-O) and 
CO(2-l) IRAM-30m spectra (centre and right, respectively) for our sample DETGs (boxcar smoothed and binned to 50 km s _ 1 ). The x-axis is the rest frequency 
for an object with expected redshift given in Table^ The grey shaded region on the spectra denotes the detected line, and the dashed lines show the baseline 
level, and ±1 a RMS level. 
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Figure 1 - continued 


4 Results 

4.1 Line detections 

We formally detect the CO(J=1-0) and CO(J=2-l) lines in 15 of 
our 17 DETGs. The CO spectra for the detections are presented 
along with SDSS 3-colour images of the sources in Figure[l] Table 


[2] lists the properties of the detected lines, and upper limits for the 
non-detections. We do not discuss the non-detections further, as the 
observations we have are not deep enough to be constraining (but 
we do include these objects in future tables for completeness). In 
all cases with detected lines the redshifts of the two CO transitions 
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Figure 1 - continued 


are consistent with one another, and with the optical spectroscopic 
redshift of the system (Table [TJ. 

The velocity width of a CO line in a galaxy depends both on 
the rotation curve of the galaxy, and the relative extent of the emit¬ 
ting components. Thus real velocity width differences between dif¬ 


ferent CO transitions may indicate that the emitting gas is concen¬ 
trated differently or that the smaller beam at higher frequencies is 
missing some gas. The velocity widths of the detected lines in our 
target DETGs are consistent in most cases, but we discuss individ¬ 
ual cases where they differ below. The velocity linewidths found for 
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Table 2. Observational parameters and derived molecular gas masses for the sample DETGs 


ID 

V S ys 

Peaki_o 

RMSi_o 

AVi_ 0 

j S„ SV L-0 

Peak 2 -i 

RMS 2-1 

AV 2-1 

J Su SV 2-1 

logio(M H2 ) 


(km s -1 ) 

(mJy) 

(mJy) 

(km s -1 ) 

(Jy km s -1 ) 

(mJy) 

(mJy) 

(km s -1 ) 

(Jy km s -1 ) 

(M 0 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

HD 

18 

11700 

32.6 

8.3 

410. 

12.3± 2.1 

68.4 

11.0 

450. 

34.2± 3.0 

9.64± 0.07 

106 

23750 

10.3 

4.3 

450. 

4.9± 1.2 

30.5 

5.7 

450. 

14.8± 1.6 

9.89± 0.11 

231 

8378 

20.4 

7.4 

358. 

7.8± 1.7 

110.0 

13.0 

364. 

41.0±3.2 

9.19± 0.10 

920 

7550 

9.5 

4.7 

391. 

3.3± 1.2 

40.0 

9.4 

286. 

12.0± 2.0 

8.65± 0.15 

967 

7936 

14.2 

4.7 

450. 

7.7± 1.3 

14.6 

5.2 

450. 

7.3± 1.4 

9.05± 0.07 

984 

10957 

17.5 

4.7 

143. 

2.5± 0.7 

31.6 

10.2 

173. 

5.l=b 1.7 

8.88± 0.13 

1004 

11675 

43.4 

10.4 

416. 

18.8± 2.6 

53.1 

17.1 

168. 

9.4± 2.8 

9.83± 0.06 

1014 

7150 

22.4 

11.0 

450. 

11.9± 2.9 

310.0 

22.0 

369. 

116.7d= 5.4 

9.21=b 0.11 

1111 

18450 

12.7 

4.2 

450. 

7.3± 1.2 

36.7 

10.2 

163. 

6.lib 1.7 

9.84± 0.07 

1232 

7650 

20.9 

7.1 

450. 

9.9± 1.9 

53.1 

19.2 

425. 

27.lib 4.5 

9.17± 0.08 

1265 

11903 

6.6 

2.7 

378. 

2.6± 0.7 

12.0 

4.8 

299. 

3.9± 1.1 

8.98± 0.11 

1267 

23419 

28.7 

6.4 

450. 

14.2± 1.8 

52.5 

14.9 

438. 

24.6± 3.6 

10.31± 0.05 

1268 

10850 

39.2 

7.8 

412. 

16.7± 2.0 

55.1 

12.9 

450. 

26.9± 3.5 

9.71± 0.05 

1372 

18530 

13.8 

2.7 

104. 

1.5± 0.4 

8.2 

6.1 

433. 

3.4± 1.6 

9.17± 0.11 

1947 

10283 

52.1 

9.1 

450. 

25.3± 2.4 

78.1 

10.0 

450. 

36.lib 2.7 

9.84± 0.04 

1030 

- 

- 

3.2 

- 

- 

- 

4.2 

- 

- 

<8.94 

1362 

- 

- 

5.7 

- 

- 

- 

- 

- 

- 

<10.19 


Notes: Column 1 lists the DETG ID number of each source. Column 2 shows the systemic velocity (in the local standard-of-rest frame) derived from fitting the 
CO(l-O) and CO(2-l) spectra. Column 3 lists the peak CO(l-O) line flux for each object, and column 4 the RMS noise reached in the observation (in «50km 
s _1 channels). Column 5 shows the CO(l-O) linewidth, and column 6 the integrated intensity of the CO(l-O) line. Columns 7-10 show the same properties as 
columns 3-6, but derived from the CO(2-l) line. Dashes indicate sources where 1mm observations were not performed. Column 11 shows the logarithm of the 
H 2 mass for each source, derived using Equation[2] For the two non detected sources Column 11 lists the 3a upper limit to the H 2 mass derived assuming a 
200km s -1 velocity width for the line. 


Table 3. Derived parameters for the sample DETGs 


ID 

(1) 

M H2 /Md 

(2) 

Mh 2 /M* 

(%) 

(3) 

M H2 /M Hi 

(4) 

Mg as /Md 

(5) 

Mgag/M* 

(%) 

(6) 

18 

172± 29 

6.2± 1.1 

- 

> 172 

>6.2 

106 

61± 15 

5.4± 1.3 

1.3ib 0.35 

107± 15. 

10.± 1.3 

231 

34±7 

4.6± 1.0 

- 

>34 

>4.6 

920 

23± 8 

0.7± 0.2 

0.9± 0.33 

50± 8. 

l.± 0.2 

967 

28±4 

0.8± 0.1 

1.9± 0.39 

44± 5. 

l.± 0.1 

984 

28± 8 

1.2± 0.4 

0.2± 0.05 

200± 8. 

9.± 0.4 

1004 

51±7 

10.2± 1.4 

0.5± 0.06 

164± 7. 

33.± 1.4 

1014 

60± 14 

4.0± 1.0 

0.2± 0.06 

315± 15. 

21.± 1.0 

1030 

<34 

<0.02 

- 

- 

- 

1111 

45± 7 

6.0± 1.0 

0.3± 0.04 

209± 7. 

28.± 1.0 

1232 

145± 27 

1.0± 0.2 

0.4± 0.09 

518± 27. 

3.± 0.2 

1265 

33± 8 

1.6± 0.4 

0.1 ±0.02 

402± 8. 

19.± 0.4 

1267 

227ib 28 

23.4± 2.9 

2.0± 0.25 

344± 29. 

35.± 2.9 

1268 

467ib 55 

7.5± 0.9 

1.8± 0.23 

730± 56. 

12.± 0.9 

1362 

<90 

<0.15 

- 

- 

- 

1372 

23±5 

3.2± 0.8 

0.1 ±0.02 

376± 6. 

52.± 0.8 

1947 

112± 10 

8.7± 0.8 

0.4± 0.04 

373± 11. 

29.± 0.8 


Notes: Column 1 lists the DETG ID number of each source. Columns 2, 3 and 4 contain the ratio of the molecular hydrogen 
mass to the dust mass, stellar mass and HI mass, respectively, while Columns 5 and 6 contain the ratio of the total gas mass 
(Hu- H 2 ) to the dust mass, and the stellar mass respectively. Dashes indicate objects where observations were not available, 
where the lower limit to the gas-to-dust and gas-to-stellar mass ratios are given. 


our sources (upto 450 km s -1 ) are consistent with the expectation 
for ETGs of this mass from the CO Tully-Fisher relation; e.g. |Davis] 
|et al.|j2011b| . 

In DETG1372 the CO(l-O) line is very much narrower than 
the CO(2-l) line. This is hard to explain physically, and noting the 
low significance of both line detections (<4cr) we suspect that one 
or both may be spurious detections. Treating the H 2 mass derived 
for this source as an upper limit in our analysis does not affect our 
conclusions, as this source is extremely HI dominated. 

DETG0920, DETG1004, DETG1014 and DETG1111 have a 


significantly smaller CO(2-l) velocity width than the CO(l-O) ve¬ 
locity width, and thus either the CO(2-l) emitting gas is more cen¬ 
trally concentrated, or the smaller beam of the telescope at 230 GHz 
causes us to miss CO(2-l) emission at higher velocities. These ob¬ 
jects do not seem to share any distinctive physical properties (e.g. 
their morphologies, starburst ages etc. seem to span the whole range 
found in the underlying population; |Shabala et al.||2012| ). [Davis 


et al. 


|2013| > studied the extent of the gas in ATLAS 3u ETGs, and 


found that the gas typically extends to half the effective radius. The 
objects with discrepant CO(2-l) velocity widths have optical effec- 
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tive radii of between 15-25", suggesting gas extents of 7.5-12.5". 
As these objects are substantially more disturbed than a typical 
ATLAS 3D object, their gas may be even more extended, making 
it feasible that an 1beam could miss some flux (see the left col¬ 
umn of Figure [T]). Beam effects are thus the more likely cause of 
this discrepancy. 

In principle, the ratio of the CO(l-O) and CO(2-l) lines can 
be used to estimate the excitation temperature of the molecular gas, 
but the significantly smaller beam size of the IRAM-30m at the 
higher frequency of CO(2-l) acts as a systematic error in such a 
measurement. Thus the line ratio we observe is affected by both 
the excitation temperature and the spatial distribution of the gas. If 
the observed CO emission was to fill the telescope beam at both 
frequencies (and the CO is not sub-thermally excited) we would 
expect a line ratio of one (measured in antenna temperature units). 
However, if the CO emission is compact compared to the beam then 
the measured intensity in the CO(2-l) line should be larger by up 
to a factor of 4 (as the beam at such frequencies covers a 4 times 
smaller area). In all but one of our objects the CO(2-1)/CO(1-0) 
antenna temperature ratio is less than 4, suggesting the total CO ex¬ 
tent is small compared to our beam. It should be noted however that 
sub-thermal excitation of CO lines would decrease the intensity of 
the higher frequency transition (e.g. |Braine & Combes|1992) , so a 
firm statement about the extent of the gas is not possible without re¬ 
solved data. DETG1014 is the exception, having a CO(2-l)/CO(l- 
0) ratio of «7.3, suggesting a high excitation temperature may be 
required in this object. 

4.2 Gas masses 
4.2.1 Hi 


The galaxies in our sample were searched for 21cm HI emission at 
the Arecibo telescope, as detailed in Allison et al., in prep, (paper 
two of this series). All but two of our objects were detected, and in 
this paper we use the HI masses (Mjj T ) as described in Allison et al. 
(see that work for full details). NGC 4797 (DETG 1232) already has 
an HI mass measurement in the ALFALFA 40% catalogue ( [Haynes | 
', which reports that Mjj T = 3.74 =b 0.54 x 10 9 M©, and 


et al. 


2011 


we use this value here. 


4.2.2 H 2 

We estimate H 2 gas masses for our CO detections in the standard 
manner, using the following equation 

Mh 2 =2m H ^-XcoD 2 L j S v SV, (1) 

where ttih is the mass of a hydrogen atom, A is the wavelength, kb 
is Boltzmann’s constant, Dl is the luminosity distance, f S v 8V is 
the integrated CO flux density and Xco is your CO-to-H 2 conver¬ 
sion factor of choice. For CO(l-O) this formula reduces to 

"»■ = 3.83 x 10-'Xco (^) a (3^1) ■ m 

Here we use a galactic Xco factor of 3x 10 2 ° cm 
km s -1 ) -1 , equivalent to aco= 4.4 M©(K km s -1 pc 2 ) -1 
|man, Snell & Schloerb|1986~| . An Xco more appropriate for low- 
metallicity gas would produce even higher H 2 masses. We return to 
discuss this assumption later, in the light of our derived gas-to-dust 
ratios (see Section [5*. 1.1) . 

We find that the detected galaxies are very gas rich, having 
H 2 masses between 4x 10 8 and 2x 10 10 M 0 . The value derived for 
each object is listed in Table [2] W e discuss the gas rich nature of 
these objects further in Section [5T| 


- 2 (K 

(Dick- 



Figure 2. A histogram of the measured gas-to-dust ratios in our sample 
DETGs, taken from Table [3] Also shown as shaded areas are the typical 
gas-to-dust ratio found in massive spiral galaxies, and early type galaxies 
(from the Herschel Reference Survey; [Smith et al.|2012), and the measured 
gas-to-dust ratio in the Large Magellanic Cloud (LMC; [Rom an-Du val et al.| 

|20T0) . 


4.3 Gas-to-dust and gas-to-stellar mass ratios 

In Table[I]we presented dust and stellar mass estimates for these ob¬ 
jects, derived from multi-wavelength SED fitting, following K13. 
In Table [3] we present the molecular gas-to-dust and stellar mass 
ratios. We also combine the molecular gas masses from this work 
with atomic gas masses from Allison et al. to tabulate molecular-to- 
atomic gas ratios and total gas-to-dust (and gas-to-stellar mass) ra¬ 
tios for our sample. Our galaxies break into two rough groups, with 
^1/3 of the sample being very molecular gas rich (i.e. Mh 2 /M^j 
> 1; DETG106, 967, 1267 and 1268) while the rest are atomic gas 
dominated. 

We find molecular-to-stellar mass ratios (Mh 2 /M*) of between 
0.7 and 24%, and total gas-to-stellar mass ratios (M gas /M*) be¬ 
tween 1 and 52%. The gas-to-dust ratios (M gas /Md) of our sample 
vary widely, between values of ^50 and 750 (see Figure [2]). 

4.4 Gas phase metallicity estimates 

The wide range of gas-to-dust ratios found in our objects imply a 
large range of gas-phase metallicities. If we assume the same frac¬ 
tion of the major condensable elements are in solid form as in the 
Milky Way (MW), we can use the results of [Draine et al.| ( |2007) to 
estimate the gas phase metallicity: 


12 + log 10 (O/H) — 12 + log 10 


f 4.57088 x 10~ 2 \ 
\ M gas /Mdust J 


(3) 


The numerical factor in the above equation is set so that 
at a gas-to-dust ratio of 100, one would predict solar metallic¬ 
ity gas 12 + log 10 (O/H) = 8.66 ([Asplund et al. 2004). Using 
the above equation we find that our objects are predicted to have 
gas phase metallicities between a tenth and twice solar metallicity 
(12 + log 10 (O/H) = 7.8 to 9.0). We discuss this result further in 
Section 15^21 


5 Discussion 

5.1 Gas masses in DETGs 

The bulge-dominated DETGs studied in this work are very gas 
rich. Each detected object has a molecular gas mass >4 x 10 s 
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M 0 (see Section 4.21, and the most gas rich galaxy has ^2xl0 10 


M 0 of molecular gas. In comparison the ATLAS 3D survey |Cap 


Ipellari et al.|2011| ), which included a complete and volume limited 
exploration of the gas content of 260 local ETGs (Young et al.|2011| 
Serra et al.|2012| ), only contained three galaxies with >10 y M 0 of 


molecular gas. 

The selection criteria we used to define our sample likely ex¬ 
plain the gas-rich nature of these objects. The galaxies show large 
dust lanes in shallow SDSS imaging and are the brightest DETGs 
in the sky in the far infrared, already suggesting they must contain 
reasonable amounts of ISM material. The selected DETGs are a rea¬ 
sonably rare population of very gas rich bulge-dominated galaxies, 
which are the ‘tip of the iceberg’ for a larger population of gas¬ 
hosting ETGs, such as that explored by ATLAS 3D . They are a good 
sample for tracing the minor merger process however, as the optical 
disturbances (see Kaviraj et al. 2012, Shabala et al. 2012 ) and large 
remaining gas masses both suggest the merger happened recently. 

5.1.1 X co 


The molecular masses used here assume a Galactic Xco factor (the 
conversion between CO flux and H 2 ). Xco has been shown to vary 
as a function of metallicity (e.g. |Wilson|1995[|Bolatto et alipOOS] 
|Leroy et al.|2011[ |Sandstrom et al.|2013) and in strong starbursts 
(e.g.|Solomon et al.|1997||Downes & Solomon|1998|i. See|Bolatto,| 
|Wolfire & Leroyfpigf for a review of this issue. 

Our objects are not starbursts, and do not have suitable strong 
emission lines that allow us to directly estimate the ISM metallicity. 
However, as shown in Section |4~4| the gas-phase metallicity may be 
sub-solar in a large number of these objects (see Section [5^2] ). Cal¬ 
ibrations between gas-phase metallicity and Xco do exist, but we 
are unable to apply these here as our proxy for metallicity (the gas- 
to-dust ratio) itself depends on Xco • We proceed from this point 
using the molecular gas masses estimated from a fixed, solar metal¬ 
licity Xco , but remain aware that we may be underestimating the 
amount of molecular gas present in these systems (by a factor of up 
to 20 in some cases; |Leroy et al.|2007| ). Total gas masses are bet¬ 
ter constrained, especially in the ^2/3 of our sample galaxies that 
are atomic gas dominated (see below). If substantial amounts of gas 
exists outside our telescope beam (which covers the inner &22" 
of these galaxies) then this would also lead to underestimating the 
molecular gas mass present. 


5.1.2 Molecular vs atomic gas 

In Section |4.3| we presented the molecular-to-atomic gas ratios in 
our objects, finding that ^1/3 of the sample are very molecular gas 
rich (i.e. Mh 2 /M Hi > 1; DETG106, 967, 1267 and 1268) while 
the rest are atomic gas dominated. The molecular-to-atomic gas ra¬ 
tio does not correlate with the gas-to-dust ratio in this sample of 
objects, or with their redshift (as would be expected if beam effects 
cause us to miss CO emitting gas in the nearest objects). In the 
ATLAS 3D survey, molecule dominated galaxies also made up ^1/3 
of the population, with a similar median molecular-to-atomic gas 
ratio (^0.3-0.5). A Mann Whittingly U test is unable find differ¬ 
ence in the distributions of molecular-to-atomic gas ratios, suggest¬ 
ing that in our DETGs the gas is in a reasonably similar physical 
state to that in the ATLAS 3D objects (despite the large total gas 
mass present). 


5.1.3 Gas to stellar mass fractions 


As discussed previously our DETGs contain large amounts of gas 
and thus, despite being very massive objects, their gas masses are 
reasonable proportions of the total baryonic mass of the system. We 
find (see Table[3} molecular-to-stellar mass ratios (Mh 2 /M*) of be¬ 


tween 0.7 and 24%, and total gas-to-stellar mass ratios (M gas /M*) 
as high as 52% in our most extreme object (DETG1372). The re¬ 
cent merger these objects have experienced causes these objects to 
be as gas rich as a typical Sb or Sc type spiral galaxy in some cases 
(e.g. [Young & Scoville|199l[|Boselli et al.|2014| . 

5.2 Dust & gas accretion in minor mergers 

As discussed in detail above, the galaxies in our sample are ex¬ 
pected to have undergone a recent minor merger. For optically red, 
massive bulge-dominated galaxies like these it is likely that the 
massive progenitor was gas poor, and the minor companion brought 
in the vast majority of the gas and dust into the system (see e.g. 
|Davis et al.|2011a| for a dynamical study showing an external origin 
for the gas in >50% of such galaxies, and | Young et al.|2 014 for an 
in-depth discussion of this hypothesis in the context of the colours 
of local ETGs). In addition, if a substantial gas mass were present 
in these objects before the merger, then the acquisition of more gas 
would likely cause a build-up of gas in the galaxy’s centre, trigger¬ 
ing a nuclear starburst. We do not see any evidence for such nuclear 
starbursts in our objects (which would be obvious in e.g. the SDSS 
spectra), suggesting once more that these objects were gas poor be¬ 
fore their recent merger. 

Further evidence for the external origin of this gas comes from 
the measured gas-to-dust ratios (listed in Table [5] and shown in 
Figure [2) and gas-phase metallicities (discussed in Section |4.4| ). 
Gas which has been produced through secular processes (e.g. stel¬ 
lar mass loss) would be expected to have a high metallicity, and 
thus a larger ratio of dust per unit gas mass than ISM acquired from 
mergers with smaller companions. Figure [2] indicates the typical 
gas-to-dust ratio found in massive spiral galaxies, and early type 
galaxies from the Herschel Reference Survey ( [Smith et al.|2012| ). 
Also shown is the measured gas-to-dust ratio in the Large Magel¬ 
lanic Cloud (LMC; [Roman-Duval et al.||2010| ), which has a sub¬ 
solar metallicity. Our sources have a wide range of gas-to-dust ra¬ 
tios, suggesting the gas was not formed in situ in the majority of 
these DETGs. 

The top panel of Figure [3] shows this in another way, plotting 
our DETG gas-to-dust ratios (and predicted gas-phase metallicities) 
as a function of stellar mass (black points). We also show the ex¬ 
pected gas-to-dust ratio as a function of galaxy mass as a black 
line (assuming the gas-to-dust vs metallicity relation from Section 
|4.4| and using the mass metallicity relation for low mass and dwarf 
galaxies; [Lee et al.|2006| ). We choose to display the dwarf galaxy 
relation here, as we will utilise it in what follows. Our extrapolation 
of this relation to high galaxy masses is not strictly valid, as above 
stellar masses of «10 10 M 0 observed mass-metallicity relations are 
found to turn over (being flatter at high stellar masses; e.g. |Tremonti| 
|et al.|2004| , but this would not change the results we derive from 
this relation, as discussed below. Mass - metallicity relations are 
notably hard to calibrate, with large systematic offsets between dif¬ 
ferent calibrations (e.g. [Kewley & Ellison||2008| ). As our galaxies 
have very similar masses, no matter which relation we take as in¬ 
put here, the spread in gas-to-dust ratio can only be explained by a 
varying metallicity. 

Three of our objects (DETG231, 920 and 967) lie on or below 
the black line in the top panel of Figure [3] having more dust per 
unit gas than expected even for their current galaxy mass. One ob¬ 
jects has no HI observations, and so its gas-to-dust ratio is likely to 
be higher, but the other two objects are genuinely dust rich. These 
galaxies have some of the most relaxed dust lanes in our sample 
(see Figure 1), and thus perhaps the gas is internally generated (and 
the tidal debris seen in deep imaging was from a gas poor minor 
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Figure 3. Top panel: The gas-to-dust ratio (plotted on the left axis) and the 
predicted gas-phase metallicity (right axis; see Section [4~4) of our sample 
DETGs plotted as a function of galaxy stellar mass (black points). Shown 
as a black line is the expected gas-to-dust ratio as a function of galaxy mass 
(using the mass metallicity relation for low mass and dwarf galaxies; [Lee] 
|et al.|2006) , with the la scatter from the mass-metallicity relation shown as 
dashed lines. Each point where the measured gas-to-dust ratio is higher than 
expected is connected to its expected minor merger partner mass by a grey 
dashed line. Bottom panel: Histogram of the predicted merger mass ratios 
inferred as in the top panel. Also shown as a red histogram is the gas-rich 
merger mass fractions predicted from the SAMs of |Lagos et al. 1(2014) . 


merger). Alternatively the gas may have been brought in a suffi¬ 
ciently long time ago that signatures of a major merger have dissi¬ 
pated, or the metallicity of gas brought in through a minor merger 
has had time to increase. From this point on we assume that the 
gas in these objects was brought in through a equal mass major 
merger, but entirely removing these objects from our analysis does 
not change our conclusions. 

If we assume, as discussed above, that the ETG present before 
the merger was completely devoid of gas, then all the gas and dust 
currently present must have been accreted from the minor merger. 
We can use the observed gas-to-dust ratio to infer the likely mass of 
the accreted smaller object, and the mass ratio of the merger. Figure 
[3] demonstrates how we do this. In the top panel the grey dashed 
lines connect the measured gas-to-dust ratio with the mass of the 
system that would be expected to have such a ratio. As many of 
the predicted stellar masses correspond to dwarf galaxy sized ob¬ 
jects, we use the mass-metallicity relation for low-mass and dwarf 
galaxies from |Lee et al.|2006| to make these predictions. The bottom 
panel of Figure [3] shows the predicted mass ratio of the merger as 


a black histogram. The median predicted mass ratio for the merger 
would be ^40:1. 

We note that the progenitor masses we derive have a high level 
of uncertainty. For instance, they are sensitive to our choice of mass- 
metallicity relation. If instead of using the result for dwarf galaxies 
(from Lee et al. 2006} we had used the relation more applicable to 
massive galaxies from |Tremonti et alT| ( |2004| ), then our derived pro¬ 
genitor masses would be lower (by up to ^0.8 dex at stellar masses 
around 10 9 M©; less for smaller satellites). The mass-metallicity 
relation has also been found to have a secondary dependance on 
star-formation rate (e.g. |Lara-Lopez et al.[|2010[ |Mannucci et al.| 
[2010) ). Our predicted progenitors are predicted to be gas rich (see 
below), and are thus likely to have been highly star-forming, and 
thus have a higher metallicity than expected for their stellar mass. 
This would lead us to underestimate the stellar mass of the progen¬ 
itor (by upto ~0.5 dex; |Mannucci et al.|2010| ). Another uncertainty 
is the state of the ETG before the merger. If we relax the assump¬ 
tion that the parent ETG was gas poor, and instead allow it to have 
some (internally produced) gas and dust which has now mixed with 
the accreted gas, then the progenitor masses we derive would go 
down. Overall, we caution that the masses derived are indicative 
only, but go on (below) to find they do agree reasonably well with 
the predictions from simulations. 

In the bottom panel of Figure [3] we show the predicted distri¬ 
bution of mass ratios for z =0 gas-rich mergers with ETGs, taken 
from the simulations of Lagos et al. ( 2014). Lagos et al. ( 2014} pre¬ 
sented a detailed study of the atomic and molecular gas contents of 
model ETGs using a semi-analytic model (SAM) of galaxy forma¬ 
tion in a ACDM universe. They extensively compared the properties 
of model ETGs with the results of ATLAS 3D . Their model repro¬ 
duces well the fraction of ETGs as a function of stellar mass in the 
local Universe, as well as the HI and H 2 content of those ETGs. 
This makes the Lagos et al. model a good tool to look into model 
ETGs that underwent a gas-rich merger. 

In order to perform a fair comparison with our observations, 
we select ETGs in the model at z — 0 (those with bulge-to-total 
stellar mass ratios>0.5) that have stellar masses > 10 10 M© and had 
a gas-rich (gas fraction of the merger system (Mjj j+MH 2 )/M s t a r > 
0.1) merger in the last ~0.4 Gyr. If we now analyse the nature of 
this sample of ETGs, we see that as expected the models predict 
these gas-rich mergers to be essentially all minor mergers (with the 
mass ratio distribution shown in the bottom panel of [3}. The other 
prediction is that most of the gas the model ETGs have after the 
minor merger was brought in by the satellite galaxy, typically in¬ 
creasing the total ISM mass of the coalesced system by a factor of 
«9, consistent with our assumption above. 

Comparing the model with our data, we find that our esti¬ 
mated companion masses follow the expected trend, and a Mann- 
Whitney U-test is unable to reject the null-hypothesis that our pro¬ 
posed mergers are drawn from the same underlying mass ratio dis¬ 
tribution as predicted from the SAMs (any differences are at less 
than a la level). The mass ratio of the gas-rich mergers is a com¬ 
pletely new constraint to the semi-analytic models, and can clearly 
reveal if the models may be having problems with the predicted gas 
contents of different galaxy populations. The fact that the Lagos et 
al. model described these observations well is encouraging as none 
of these observations were used to construct the model. The Lagos 
et al model, despite its simplifications, does seem to capture at least 
some of the behaviour seen in real galaxies. 

In Figure[4]we plot the gas fractions the gas-rich dwarf galaxy 
is predicted to have had before the merger as a black filled his¬ 
togram. The gas fractions of some of the accreted satellites are high, 
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Figure 4. Histogram of the predicted gas fractions in the minor merger com¬ 
panion, assuming all the gas mass in our DETGs came from a minor merger 
with a dwarf companion (whose mass is derived as in Figure[3j. Also shown 
as a shaded histogram are the observed gas fractions in dwarf galaxies from 
jHuang et aL](2012) , and as a red histogram is the gas fractions predicted for 
gas rich mergers onto ETGs in the SAMs of |Lagos et al. 1(2014) . 
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Figure 5. Histogram of the star formation efficiency (SFE=SFR/Mh 2 ) 
in our DETGs (black histogram) and ATFAS 3D ETGs (grey shaded his¬ 
togram). The average SFE of the spiral galaxies from |Kennicutt| (1998) is 
also shown as a black vertical line, with a grey bar showing the standard 
deviation of their SFEs around the mean. 


with the median accreted galaxy having ~10 times more gas mass 
than stellar mass. The error bars on this estimate are large, with a 
median error bar being displayed on the bottom right (we include 
the errors in gas and dust masses, and propagate these through, in¬ 
cluding a source of error from the scatter of the mass-metallicity re¬ 
lation). In Figure [4] we also plot the observed distribution of dwarf 
galaxy gas mass fractions from |Huang et al.| ( |2012| ), shown as a 
grey shaded histogram. We also show as a red histogram the pre¬ 
dicted gas fractions of z =0 gas-rich mergers onto ETGs, again taken 
from the semi-analytic models (SAMs) of |Lagos et al.H2014| ). The 
predicted and observed gas fractions (from |Huang et al.[2012| and 
Lagos et al. 2014 respectively) are found to agree well with one 
another. 

As Figure [4] shows, our predicted distribution of minor com¬ 
panion gas mass fractions are skewed to larger values than found in 
the sample of local dwarf galaxies. This can likely be explained by 
our strong selection function towards detecting more gas-rich ob¬ 
jects, ensuring we are only selecting the most gas rich minor merg¬ 
ers in the local universe. It should be noted that the observed dwarf 
galaxy gas fractions actually only include atomic gas, thus the total 
gas mass is somewhat higher, decreasing any discrepancy. Dwarf 
galaxies as gas-rich as those required here do exist in the study of 
|Huang et~ak] ( |2012| ), but are rare («7% of the observed dwarfs have 
gas to stellar mass fractions greater than 20). They also found «2% 
of dwarfs had gas masses greater than 100 times their stellar mass. 

5.3 Suppression of star-formation in minor mergers 

In Table [l] we presented the star-formation rates of the sample 
DETGs, estimated via SED fitting (using the MAGPHYS code). The 
derived values vary between 0.01 and 1.3 M 0 yr -1 . These are rea¬ 
sonable star-formation rates for ETGs, well within the range found 
within the ATLAS 3D sample (0.008 - 18 M 0 yr _1 ; D14). 

|Davis et~ak] < |2014[ hereafter D14) reported that the star- 
formation was suppressed in ATLAS 3D ETGs (compared to late- 
type galaxies), with their sample galaxies lying on average a fac¬ 
tor of ~2.5 below the Kennicutt-Schmidt relation (e.g. Kenni- 
|cutt|1998] ). We are unable to place our galaxies on the Kennicutt- 
Schmidt relation, as we have no estimate of the size of the star¬ 
forming region in these objects. We can however estimate the star- 


formation efficiency (SFE; defined as star formation per unit gas 
mass), using our gas masses and derived star-formation rates. In 
what follows we define the SFE as star-formation rate per unit 
molecular gas mass (SFR/Mh 2 ) rather than total gas, because 
molecular gas is the phase which is most closely linked to star for¬ 
mation (e.g. |Bigiel et al.|2008| >. Including Hi (or using a variable 
Xco) in our SFEs would make the SFEs even smaller (or equiva¬ 
lently, make depletion times longer). 

Figure [5] shows the SFE in our DETGs (black histogram), 
and the ATLAS 3D ETGs (grey histogram). Also shown (as a black 
line/grey bar) are the spiral galaxies of |Kennicutt| ( |1998] ), which 
have been found to have average SFEs of «1.5x 10 y yr 1 (see 
also e.g. |Bigiel et al.|20lT] ). The distribution of SFE in our DETGs 
overlap with the distribution found in ATLAS 3D , but the medians 
are substantially different. Our DETG sample objects have median 
SFEs of«9 x 10 11 yr 1 , while the median for the entire sample of 
ATLAS 3D ETGs is closer to ^4 x 10“ 10 yr -1 . A Mann-Whitney U 
test finds that these two SFE distributions are unlikely to be drawn 
from the same parent distribution at a confidence level ^>3 a. Thus, 
if star formation were to continue on in a steady state at the same 
level as today, the gas reservoir in our objects would take a me¬ 
dian time of wll Gyr to be depleted. We note that the median SFE 
in our sample would be even lower if we use a metallicity-derived 
Xco conversation factor. 

We checked that our results were not dependant on the method 
used to estimate star-formation rates by also calculating monochro¬ 
matic star-formation rates from Wide-field Infrared Survey Explorer 
(WISE; [Wright et al.|2010| ) 22 pm emission alone, and the 22 pm 
data combined with observations of far-ultraviolet emission from 
the Galaxy Evolution Explorer (GALEX; |Morrissey et al. 


We used the same procedure as described in |Davis et ah 


2007). 


2014), 


taking 22 gm aperture magnitudes and foreground extinction cor¬ 
rected FUV magnitudes for the DETGs from the data compiled in 
K13. We note that the aperture values used include the entire galax¬ 
ies, allowing us to properly subtract the circumstellar emission from 
old stars, which can be important at 22 pm (see |Davis et al.|p014) 
for a full discussion of this process). These three estimates of the 
star-formation rate agree well, apart from at low SFRs where the 
22pm indicator gives slightly higher SFRs (at a la level). What- 
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ever the method used to estimate the SFR, we find our DETGs have 
very low SFEs, with median depletion times in excess of 6.6 Gyr. 

This low star-formation efficiency in remnants of recent gas- 
rich minor mergers is puzzling. Major mergers are usually associ¬ 
ated with strong starbursts at low redshift (e.g. Arp 220), with ex¬ 
treme star-formation rates , an d fast d epletion times (some as low 
as 10 7 yr; e.g. Gao & Solomon|2 004). This is usually explained as 
a dissipative collapse of the gas in the galaxy to the centre, driv¬ 
ing gas densities up, which then results in a nuclear starburst (e.g. 
[Mihos & Hernquist|1996| ). The three objects in our sample which 
may have had major mergers (discussed in Section |5^2] do lie at the 
higher end of our SFE distribution, but are not strong outliers. Mi¬ 
nor mergers are less violent events, but these too might be expected 
to increase the SFE by some factor (see e.g. |Saintonge et al.|20T2] 
who report evidence for such a relation with mergers of all types). 
The DETGs studied here do not have strong AGN, supporting the 
idea that there is a large time delay (in merging systems) between 
the onset of star formation and AGN activity ( [Shabala et al.|2012[ 
Kaviraj et al. 2014). We thus consider that AGN feedback is un¬ 
likely to be important in suppressing the star formation in these 
objects. 

Studies such as Saintonge et al. ( 20l2]> and[Kaviraj ( 2014a]) do 
find that the presence of a large bulge reduces the SF E (a conclu- 
sion supported by the low SFE of ATLAS 3D ETGs in Davis et al. 


|2014| ), and it has been posited that this effect may arise because the 
presence of a bulge keeps gas more stable against star formation 
(e.g. Martig et al. 2009 2013). In our sample, the galaxies with a 
more ‘normal’ SFE do tend to have strong disks, but so do some of 
the very suppressed objects. It is also unclear why this sample of 
recent minor merger remnants would show even greater SFE sup¬ 
pression than normal ETGs (both those from ATLAS 3D , and the 
massive galaxies in the |Saintonge et al.|2012| sample), as they likely 
have similar bulge masses. A relation between stellar mass and star- 
formation efficiency has been shown by |Saintonge et al.| |2012| ); 
Bothwell et al. (2014], but our objects lie well off this relation. Sta¬ 
bility enhancement by the bulge would also likely be most efficient 
on regular gas reservoirs which have been in place for some time. 

One possible explanation is that we have caught our objects 
in a special phase of their evolution, where gas is still streaming 
towards the centre (potentially being stabilised by these inflowing 
motions, e.g. |Meidt etaL||2013[ or additional shocks/turbulence). 
The size of the star-forming region in this type of galaxy has been 
observed to decrease as the merger proceeds (Shabala et al., in 
prep.), supporting the idea the gas may be inflowing. Once gas 
densities in the centre of these objects increase a starburst may be 
triggered, which would lead us to no longer classify these galaxies 
as DETGs. Support for such a picture comes from observations of 
star-formation in tidal tails produced by minor mergers and ram- 
pressure stripping, where the SFE is also found to be low (e.g. 
Knie rman et al . 2013, Jachym et al. 2014). In addition, observa¬ 
tions of interacting systems/compact groups suggest that if the ISM 
is kept warm by shocks and accretion induced turbulence, this may 
result in suppressed star-formation efficiencies (e.g. [Alatalo et al.| 
|20l4| | Appleton et al.|2014| ). 

Conversely, many of the CO line profiles shown in Figure[T]are 
(within our ability to measure them) symmetric and show double 
horns (that normally result from gas rotating regularly in the galaxy 
potential). Paper II presents more in depth analysis of the HI line 
shape, which may be more sensitive to dynamical perturbations (be¬ 
cause dynamical times are longer in large HI disks). The dust disks 
visible in the optical images of these galaxies for the most part also 
seem regular. These two pieces of evidence argue against the inter¬ 


pretation above, where the ongoing merger is strongly dynamically 
perturbing the gas reservoir, but without resolved gas imaging it is 
hard to conclude what the correct physical explanation for the low 
SFE of our DETGs may be. 

6 Conclusions 

In this work we used the IRAM-30m telescope to detect CO 
emission in a sample of bulge-dominated galaxies which have large 
dust lanes visible in optical imaging, and far-infrared fluxes avail¬ 
able from the Herschel space telescope. These systems show distur¬ 
bances in their optical light distribution that suggest they have had 
minor mergers in the recent past, which have supplied their gas and 
dust. 

We firmly detect 15 of our target objects in at least one CO 
line, and determine that these galaxies are very gas rich, having Fb 
masses between 4x 10 8 and 2x 10 10 Mq. We use these molecular 
gas masses, combined with atomic gas masses (from Paper II in this 
series) to calculate gas-to-dust and gas-to-stellar mass ratios (using 
dust and stellar masses derived from SED fitting following K13). 
We find molecular to stellar mass ratios (Mh 2 /M*) of between 0.8 
and 23%, and total gas to stellar mass ratios (M gas /M*) between 1 
and 51%. The gas-to-dust ratios (M gas /Md) of our sample objects 
vary widely, between values of ^50 and 730 (see Figure [2]. We 
discuss our assumption of a fixed Xco for these objects given this 
variation, and conclude that we may be underestimating the molec¬ 
ular gas mass in many of our objects. Any such underestimation 
only serves to strengthen the results that follow. 

The large range in gas-to-dust ratios in our objects implies a 
wide range of gas-phase metallicities, much larger than predicted 
by the scatter in the mass-metallicity relation. This suggests that 
the gas in these objects has indeed been accreted through a recent 
merger with a lower mass companion. We use the gas-to-dust to 
metallicity, and the mass-metallicity relations to calculate the im¬ 
plied minor companion masses and gas fractions, assuming the gas 
and dust have been accreted in their current proportions. We find the 
likely dwarf companions had masses as low as «10 7 M 0 and that 
the implied merger mass ratios are consistent with the expectation 
from simulations. The median predicted mass ratio for the mergers 
is ^40:1. The accreted satellites are predicted to be very gas rich, 
but are again consistent with being drawn from the population of 
known dwarf galaxies. 

In Section [53] we showed that (no matter which SFR indicator 
is used) our sample objects have very low star-formation efficien¬ 
cies, taking on average 6.5 Gyr to deplete their gas reservoirs. This 
cannot be purely a result of morphological quenching mechanisms, 
as this level of suppression is greater than found in the ATLAS 3D 
sample of ETGs. We discuss mechanisms that could cause such a 
suppression, include dynamical effects caused by the minor merger, 
but with the current data available are unable to reach a firm con¬ 
clusion. 

These results clearly motivate further study of the minor 
merger process, and its effect on star formation in the remnant. 
Resolved observations will be required to understand the config¬ 
uration of the gas in these objects and enable resolved studies of 
their star-formation efficiency. Simulations of gas-rich minor com¬ 
panions merging with gas-poor ETGs will help determine if the ob¬ 
served suppression of star formation is caused by the merger itself, 
and thus if it is a more general process that has importance across 
the Hubble sequence. 
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